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Abstract: We update the values of the r\-rj mixing angle and of the rj' gluonium content 
by fitting our measurement Rj, = BR(</> — > 7/7)/BR(</> — > 777) together with several vector 
meson radiative decays to pseudoscalars (V — > P7), pseudoscalar mesons radiative decays 
to vectors (P — > V7) and the rf — ► 77, 7r° — > 77 widths. From the fit we extract a gluonium 
fraction of Zq = 0.12±0.04, the pseudoscalar mixing angle ipp = (40.4±0.6)° and the <f>— to 
mixing angle ipv = (3.32 db 0.09)°. Zq and ipp are fairly consistent with those previously 
published. We also evaluate the impact on the rj gluonium content determination of future 
experimental improvements of the rf branching ratios and decay width. 
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1. Introduction 

The rj meson, being almost a pure SU(3)fl avour singlet, is considered a good candidate 
to host a gluon condensate. The question of a gluonium component in the rj meson has 
been extensively investigated in the past but it is still without a definitive conclusion 
Q. We extract the rj gluonium content and the rj-rj mixing angle in the constituent 
quark model according to the Rosner |^] approach with the modifications introduced in 
ref. || as described in the following. We use the same method of ref. p]]; in addition, we 
also introduce in the fit the tt° — > 77 and rj — ► 77 branching fractions according to the 
prescription of ref. ||]. This method relates our measurement of the ratio <p ~^ v'l arm 
<f> —¥ ryy branching ratio (BR), = BR(0 — >• r]'j)/BH((f) — > 777) Q, to the rj' gluonium 
content and to the rj, rj' mixing angle. The same quantities were extracted in our previous 
analysis || with some assumptions. This has given rise to some objections from refs. Q 
and 0. Here we give an answer to these objections and we repeat the fit taking into account 
their comments. Then we repeat the fit with recently updated experimental results. 

The rj and rj' states can be represented in the base \N) = (\uu) + \dd^)/\^2, \S) = \ss) 
and \G) = | gluonium) as: 

\rj') = cosipc simpp \N) + cosipc cosipp \S) + sin-^G \G) (I- 1 ) 
\rj) = cos i/) P \N) - smip P \S) (1.2) 

where tpp is the rj-rj mixing angle and Zq = sin 2 ^ is the gluonium fraction in the 7/ 
meson. According to ref. || the state \G) could be the 7/(1405): a pseudoscalar glue ball 
candidate. 
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The ratio = BR(0 — ► 7/7) /BR^ — > 777) is related to the ipp and parameters by 
the formula Q : 

i?^ = cot>pcos> G 1-— — [ — ) (1-3) 

where py and jj.^ are the momenta of the 77' and 77 meson respectively in the <f> reference 
frame, m s /fh = 2m s /{m u + md) is the constituent quark masses ratio and ipy is the 4>-uj 
mixing angle. Following ref. Q we define the constant C q = (qq p \qq ri ) as the overlap 
between the spatial wave functions of the quark-antiquark pair in the p and the 77 meson. 
Isospin symmetry is assumed exact, so that m u = = fh and the following further 
relations follow: 

C q = (qq v \ qq~u) = {qq v \ qq P ) , C s = (ss v \ ss^} , CV = (qq n \ qq u ) = (qq w \ qq p ) 

where we indicate with \qq v ) and \qq w ) the qq spatial wave function in the 77 and u mesons, 
and with \ss„) and \ss^) the ss spatial wave function in the 77 and (ft mesons. The parameters 
Z q and Z s are the ratios: Z q = C q /C n and Z s = C s /C n . In this model SU(3)fl avoU r 
breaking effects are accounted for by the different values of the effective quark masses, 
m s > m u = rrid = fh, and by Z q ^ Z s . 

In our previous analysis || the parameters Z s , Z q , ifty and m s /fh were taken from ref. 
II where BR(<^> — > 7/7) and BR(0 — > 777) were fitted together with other V — ► P7 decay 
rates (V indicates the vector mesons p, u, (ft and P the pseudoscalars 7r°, 77, 77') assuming no 
77' gluonium content. We fitted Q our measurement 

R<t> = BR(^^7) = (4 ' 77 ± °'° 9stat - ± °' 19SS/S ° X 10 ~ 3 

together with the available data || on r(77' — > 77)/r(-7r° — > 77), r(r/ — ► pj)/T(u) —* 7r°7) 
and r(77' — > u^f)/T{u — > 7r°7). The dependence of these ratios from the mixing angle ?/p 
and the gluonium content ^3 is given by the following equations: 

X^/ = sim/p cost/y?, 5^' = cosipp cosiftc 

5^- cosipa sini/p + V2^- cosipc cosipp j (1-4) 
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where f n is the pion decay constant and f q and / s are the decay constants of the isospin 
singlet states (mainly 77,77' mesons) at no anomaly limit |5|. The fit result was iftp = 
(39.7 ± 0.7)° and Z 2 G = sin^G = 0.14 ± 0.04, P{ X 2 ) = 49%. Imposing ift G = the x 2 
probability of the fit decreased to 1%. 
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In ref. a procedure similar to Q has been adopted, but they fitted also the 

gluonium component in the rj wave function that was previously fixed at zero. The result 
Zq = 0.04 ±0.09 deviates la from our value but with a larger error. In refs.J|] and j?| this 
difference was attributed to the use in our fit of the parameters Z s and Z q obtained in ref. 
Q assuming no gluonium content. However further tests of the fit procedure showed that 
Zq and ifip are marginally sensitive to large variations of Z q and Z s Here we repeat 
the fit with a larger number of free parameters, including Z s and Z q . 



2. Fit description 

In order to enlarge the set of parameters used in the fit we add to the equations ( |1.4| - |L6| ) 
the following further relations that can be derived directly from 0: 
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Notice that, differently from where the VP7 couplings are fitted, we fit directly ratios 
of partial decay widths. This allows both to reduce the parameters involved in the fit 
(two of them cancel out in the ratios), and to simplify the error treatment using quantities 
proportional (directly or inversely) to the experimental measurements. As an example, a 
ratio of T's is written as: 

T( V ' -> pry) _ BR(r/ - pry) T v , 



In this way the correlation matrix among the rf branching ratios and the decay widths can 
be used directly . The fit is performed minimising the x 2 function: 



X 



where j/j are the experimental measurements of the ratios on the left side of equations ( 1.4 - 
q), yf 1 are the theoretical predictions on the right side of the same equations 



|1.6D and fly 

and V -1 is the inverse of the covariance matrix. This last quantity is obtained summing 
the contribution from the experimental error on branching ratios, decay widths and their 
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Table 1: Correlation matrix of the rj branching ratios from PDG-2006 fit ||. 
correlations and the uncertainty from theoretical inputs. The method is the same of our 



previous fit @, but only the parameters f q /f n = 1.00 ± 0.01 and f s /f n = 1.4 ± 0.014 ]TT] 
are taken as input, while in the past also Z s , Z q , ipy and m s /fn were fixed. The parameters 
fq/fn- an d fs/fn involve only the ratio T{rj — > 77)/r(-7r° — » 77). The contribution from 
theoretical error is evaluated by standard error propagation: 

V th = A C A T 

where C is the covariance matrix of the uncorrelated parameters x\ = f q /f^ and X2 = fs/fn 
and A is: 

du th 

(A) « = ~k 

The covariance matrix V is indeed: 

V = V exp + V th 

where V exp is the covariance matrix of the data used in the fit. Particularly relevant is the 
correlation between the rj branching fractions and the decay widths shown in tab. [l|. 

The rj width and the rj — ► 77 branching ratio are 88% correlated. This is because 
the rj width is evaluated dividing the T(rj' — > 77), obtained by measuring the cross section 
a(e + e~ — > rj'e + e~), by the 7/ — ► 77 branching ratio. The fit results and the correlation 
matrix are shown in tab. |2[ The gluonium fraction Zq = 0.105 ± 0.037 is 2.8a from zero. 
In order to identify the measurements which require the presence of the gluonium in the rj 
we have repeated the fit fixing Zq at zero. The results of the fit are shown in tab. ^. The 
X 2 probability is now quite low, reflecting the 2.8<7 effect seen in the previous fit, while the 
pseudoscalar mixing angle is quite stable. 

In fig. [j] we show for each measurement the pulls of the two fits, defined as pi = 
(Vi ~ yf l )/ a Vi- The measurement which does not fit in the no- gluonium hypothesis is the 
ratio V{rj — ► 77)/T(7r° — ► 77), whose pull is less than —3, bringing the x 2 probability to 
1.1%. It becomes ~— 1 when gluonium is allowed. We have then repeated the fit without 
using the T(rj' — * j^)/T(tt° — > 77) information. The result is compared with ref. 0] in 
tab. |3[ The results of the two fits are in agreement, making it evident that the origin of 
the discrepancy with ref. Q is due to their neglecting the T(n' — > 77)/T(7r° — ► 77) datum. 
In ref. the couplings among the vectors and the pseudoscalar mesons are used in place 
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Table 2: Fit results using the PDG-2006 data || and their correlation matrix. 
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Figure 1: Pulls of the fit shown in tab. left: Zq free, right: Zq = (fixed). 

of the width ratio. The couplings are related to the partial decay width by the following 
formulae: 

In order to make a full comparison between the two methods we have performed the fit 
also using the couplings and we have obtain the same results [O ]. 
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Table 3: Comparison among the fit results without the rj — > 77/71" — > 77 measurements and the 
results of ref. [g. PDG-2006 data §] have been used in both fits. 
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Table 4: Fit results using the PDG-2008 data. 



3. Update with the recent PDG results 



In the Review of Particle Physics |13| new measurements of the p, u, rj and rj mesons have 
been included, which change slightly the partial decay widths used in the fit. Therefore 
we repeat the fit using these updated values together with our measurement. All the 
correlation coefficients among the measurements are taken into account in the fit. The 
results of the fit are shown in tab. |I] and the correlation matrix in tab. S the pulls of the 
fit are shown in fig. ^. 

The results in tab. || show that 
the gluonium hypothesis is still hi- 
ghly favoured with respect to the 
null gluonium hypothesis. Never- 
theless the fit probability is quite 
low also in the gluonium hypothe- 
sis: it goes from 17% using PDG06 
data to 5% using PDG08 data. The Table 5: Correlation matrix from the fit shown in tab. 

reason of the worsening of the fit is 

found comparing the pulls of the new fit (fig. |2|) with the previous one (fig. |l]). In partic- 
ular, the pull of the ratio T(u> — ► 777) /T(u —* 7r°7) goes from -0.93 using PDG-2006 data 
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Figure 2: Pulls of the fit using PDG-2008 data, left: Z G free, right: Z G = (fixed). 

to -1.5 using PDG-2008 while the pull of the ratio T(p — ► 777) /r(w — > 7r°7) goes from -0.14 
to +0.39. This happens because the PDG estimate of the BR(w — * 777) has changed from 
(4.9 ± 0.5) x 10~ 4 to (4.6 ± 0.4) x 10 -4 , lowering the T's ratio and worsening the pull 1 . 
The BR(u> — ► 7/7) and BR(p — > 7/7) PDG values are dominated by the measurement of 



the e + e~ — * 7/7 cross section by SND [14] as a function of in the p, u, <p mass range. From 
the measured cross section they extract the p, uj partial decay widths assuming the Vector 
Meson Dominance model and a parametrisation for the p' resonance. Some correlation is 
therefore expected between the p and the u partial decay widths which are not discussed 
in ref. fi"4|| , moreover the decay widths are model dependent. The average value reported 
by PDG-2008, BR(lj — » 777) = (6.3 ± 1.3) x 10~ 4 , is dominated by a model independent 
measurement |ff5[ and is 1.2cr away from the PDG fit. Using this value for BR(w — ► 7/7) we 



obtain a much better x probability: P(x ) = 28% in the gluonium hypothesis and 1.] 



fixing the gluonium at zero [T^J . Both gluonium content and pseudoscalar mixing angle are 
unchanged (Zq = 0.11 db 0.04, ifip = (40.6 ± 0.7)° in the gluonium hypothesis). Therefore 
we will use the average value for BR(w — ► 777) in the following. 

4. Update with the new KLOE measurement of BR(w — > 7r°7) 



The relations ( ]1.5| - [1.6| ) and (2.1-2^) are dependent from the lo — > 7r°7 decay rate. Recently 
we have improved the measurement of this branching fraction BR(w — ► 7r°7) = (8.09 ± 



0.14)% This value is about 3cr different from the PDG 2008 value: BR(w -> ir°j) = 
(8.92 ± 0.24)%. We then performed the fit using our measurement of BR(cj — ► 7r°7). 

Moreover f q /f n and f s /f w have been fixed according to ref. ||. In the exact isospin 

1 A fit without the F(r]' — > 77)/r(7r° — > 77) ratio has been performed in order to check its effect on the 
w ~~ * VI an d p — » J?7 pulls. In this case the x 2 / n df of the fit is high also in the null gluonium hypothesis, 
nevertheless the Y(u — * 777)/T(w — > 7r°7) pull is -0.68 while the F(p — > iyy)/r(ijj — > vr°7) pull is +1.1. 
In other words the BR(u — » 777) fits better while the BR(p — » 777) fits worse. The difference of the two 
measurements from the best fit is still ~ 2a, therefore the poor \ 2 with PDG-2008 data is not due to 
77' — > 77 but to the inconsistency between T(lu — > 777) and T(p — ► 777) measurements. 
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symmetry approximation the relations 
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hold, where f w and fx are the ir and K decay constants. Therefore f s /f-K = JVjJff^l. 
Using fx/ fn from lattice calculation [19] we get: 

f f 
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The results of the fit are shown in tab. || and the correlation matrix in tabQ 
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Table 6: Fit results using PDG-2008 inputs, BR(cj — > 777) from PDG direct measurement average 
and the KLOE BR(w — > ir°j) and R^. The equations (4.1) have been used for the f q / / w and f s / fir 
parameters. 

The 77-77' mixing angle and the rj' gluonium content are not substantially modified, but 
the x 2 probability is improved with respect to the previous fits. The <fi — u mixing angle, 
tpv > is also slightly changed from (3.17 ± 0.09)° to (3.32 ± 0.09)° by our new measurement 
of BR(cj — > 7r°7), being tpy directly related to the ratio T(4> — > tt°^)/T(uj — > 7r°7) (see eq. 
U). 

In fig. |H we show the 68% CL contour of the 7/ related measurements in the Z G ,ipp 
plane where the contribution of each measurement to the fit results can be appreciated. 
The best fit values of the width ratios are shown in tab. |8| together with their correlation 
coefficients in tab. U. 



5. Prospects with improved measurements 

We have shown that within the precision of the available data the crucial measurement 
sensitive to the v( gluonium content is T(r/ — > j^)/T(tt° — > 77). The theoretical framework 
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Table 7: Correlation matrix of the fit shown in tab. 6. 
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Figure 3: 68% confidence level regions of the shown measurements in the Z%,t/jp plane. 
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0.06 



Table 8: Fitted values of the T ratios. 
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Table 9: Correlation coefficients of the best-fitted values. 



used to describe 7/ — > 77 and V — > P7 transition is different. In fact while in the first 
case we have a quark-antiquark annihilation into two photons, in the second case we have 
a transition among two meson with a photon emission via the spin flip of one of the two 
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(40.5 ± 0.6)° 


(40.5 ±0.6)° 
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0.93 ± 0.02 


0.93 ± 0.03 


Z s 


0.83 ± 0.05 


0.83 ± 0.05 


ipv 


(3.32 ± 0.08)° 


(3.32 ±0.09)° 


m s /fh 


1.24 ± 0.07 


1.24 ± 0.07 



Table 10: Fit values assuming 1% error on the rj branching fractions: left: using rj — > 77/71"' 
77, right: without using rj — > 77/71° — + 77. 



quarks. Therefore it is important both to reach a sensitivity to the gluonium independently 
from the rj — > 77 decay and to measure again the 7/ — > 77 branching ratio. 

The rj' — * 77 branching ratio affects directly the rj' decay width while the rj' — * ix^ix^rj 
and the rj — > ir°ir°rj branching ratios affect the systematic errors on the measurement 
H. The increase of the precision on all main rj' branching fractions at 1% level could 
lead the gluonium contribution to ~ 4cx (see tab. |l0| ). The results shown in the table 
are obtained assuming as central value of all measurements the width ratios shown in 
tab. |8| and assigning to all rj' branching ratios a 1% error (the error on r„/ and 
has been recomputed accordingly); the correlation matrix is taken equal to the present 
measurements. The fit gives a large rj gluonium component even without using the rj — > 
77/71° — ► 77 measurement (see tab. 10, right column). 



In fig. U the expected 68% CL contour in the plane i/jp,Zq is shown. Notice that 
the rj' — ► 77 measurement would be not needed anymore to have a significative gluonium 
component. 

A big improvement in the gluonium content determination can be obtained by mea- 
suring the rj' decay width directly through the measurement of cr(e + e~ — ► e + e~7*7* — > 
e + e~rj'). The measurement of this cross section at 1% level and the 7/ — > 77 branching 
ratio at the same level of accuracy would allow to determine the rj' width: 

Y(rj -» 77) 



BR(r/' -» 77) 



with a fractional error of ~ 1.4%. In fig. |], right the 68% CL contours in the (iPp,Zq) 
parameter plane are shown for all rj related measurements. The improvement in the fit 
accuracy is evident. The gluonium contribution would be measured at ~ 5er's. 

6. Conclusions 

The origin of the difference between the KLOE result || on the rj gluonium content 
and ref. 0] is the use of the Y(rj — > 77)/r(7r° — * 77) measurement. A global fit to 
all measured V — * -P7 and P — > V7 transitions of light mesons has been performed 
extracting all the relevant parameters. The new results are slightly different from our 
previous results but confirm the presence of a significative gluonium contribution in the rj 
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Figure 4: Future prospects: one sigma confidence regions in the (ipp, Zq) plane assuming present 
fitted measurements but with reduced errors: left: with 1% error on rf branching ratios, right: as 
before plus 1% precision on T(rf — ► 77). 



meson. The origin of this contribution has been investigated and found to stem mainly 
from the T(r)' — > 77)/r(7r°77) measurement. The fit has been updated with all recent 
measurements from PDG [|l3| and some disagreement with the PDG fitted value of the 
BR(iij — ► 777) has been found. The average of the u — > 777 branching ratio is therefore 
preferred to the PDG fit and a better fit is obtained while the gluonium content and 
the mixing angle are unaffected. The result has been updated using the recent KLOE 
measurement of the to — > 7r°7 branching ratio and and the lattice results for f g /f n and 
fs/ fir assuming exact isospin symmetry. The mixing angle and the gluonium content are 
again marginally affected by this measurement while a different <p — u) mixing angle is 
obtained. 

Finally, we have estimated the impact of future improvements in the measurements of 
the 7/ branching ratios and rf width in the gluonium content determination. A 5cr evidence 
could be obtained with a 1% precision on each measurement. 
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